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Abstract 

A possible model of twin high-frequency QPOs (HF QPOs) of microquasars is 
examined. The disk is assumed to have global magnetic fields and to be deformed 
with a two-armed pattern. In this deformed disk, set of a two-armed (m = 2) vertical 
p-mode oscillation and an axisymmetric (m = 0) g-mode oscillation are considered. 
They resonaltly interact through the disk deformation when their frequencies are the 
same. This resonant interaction amplifies the set of the above oscillations in the case 
where these two oscillations have wave energies of oposite signs. These oscillations are 
assumed to be excited most efficiently in the case where the radial group velocities 
of these two waves vanish at the same place. The above set of oscillations is not 
unique, depending on the node number, n, of oscillations in the vertical direction. 
We consider that the basic two sets of oscillations correspond to the twin QPOs. 
The frequencies of these oscillations depend on disk parameters such as strength of 
magnetic fields. For observational mass ranges of GRS 1915+105, GRO J1655-40, 
XTE J1550-564, and H1743-322, spins of these sources are estimated. High spins of 
these sources can be described if the disks have weak poloidal magnetic fields as well 
as toroidal magnetic fields of moderate strength. In this model the 3 : 2 frequency 
ratio of high-frequency QPOs is not related to their excitation, but occurs by chance. 

Key words: accretion, accrection disks — black holes — quasi-periodic oscilla- 
tions — X-rays; stars 



1. Introduction 

Many quasi-periodic oscillations (QPOs) have been observed in low-mass X-ray binaries 
(LMXBs). Among them, high frequency QPOs are particularly interesting, since they will 
be attributed to the innermost regions of relativistic accretion disks and thus clarification of 
their origins will give important informations about disk structure in strong gravitational fields 
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as well as about mass and spin of the central sources. High frequency QPOs in LMXBs are 
classified into two classes: High frequency QPOs (HF QPOs) in black-hole LMXBs and kHz 
QPOs in neutron-star LMXBs (van der Klis 2004). HF QPOs of black-hole LMXBs have no 
time variation and in some sources they are observed in pairs whose frequency ratio is close to 
3 : 2. KHz QPOs in neutron-star LMXBs are also observed in some sources in pairs, but their 
frequency ratio is not 3:2, but changes with correlated time variations. It is unclear whether 
both types of QPOs come from a common origin. 

So far, many models have been proposed to describe HF QPOs of black-hole LMXBs, 
especially focusing on the 3 : 2 frequency ratio. Typical models will be the relativistic precession 
model proposed in a series of papers by Stella & Vietri (1998) and Morisink & Stella (1999), and 
the epicyclic resonant model proposed by Kluzniak & Abramowicz (2001) and Abramowicz & 
Kluzniak (2001) and extensively studied by their collaborators (e.g., Abramowicz et al. 2003a, 
b; Bursa et al. 2004; Rebusco 2004, Kluzniak et al. 2004, Horak 2008, Stuchlik et al. 2008a, 
b, Horak et al 2009 and others). In a context different from the epicyclic resonant model, 
importance of a resonance on wave excitation has been emphasized, and the warp resonant 
model has been proposed (Kato 2003, 2004, 2008, Kato & Fukue 2006, Ferreira & Ogilvie 2008, 
Oktariani et al. 2010). This model is based on resonant excitation of a set of positive-energy 
and negative-energy oscillations by interaction through disk deformation (Kato et al. 2011, 
Kato 2011b). In addition to the above models, many discoseismic modes have been studied in 
order to examine their applicability to QPOs since Kato & Fukue (1980) (see Wagoner 1999, 
Kato 2001 and Kato et al. 2008 for reviews, and Wagoner 2012 for recent work). 

In spite of many efforts by the above-mentioned work and by other reseraches, the twin 
QPOs observed in microquasars are not consistently described by one unique model (Torok et 
al. 2011). This difficulty becomes prominent by recent findings that microquasars with twin 
QPOs have extremely high spins (McClintock et al. 2011). 

In the case of kHz QPOs in neutron-star LMXBs, frequencies of twin QPOs and their 
time variation seem to be well described by assuming that kHz QPOs are two-armed vertical 
p-mode oscillations trapped in the innermost region of disks (Kato 2011a,c,; 2012a,b,c). In the 
case of twin HF QPOs of black-hole LMXBs, however, a direct application of trapped vertical 
p-mode oscillations seems not to well describe observations. 

One of important clues to explore the origins of the twin HF QPOs of black-hole LMXBs 
will be why they appear only in or near to the very high state (steep power-law state) (Remillard 
2005), where the disk consists of both a power-law component and a thermal disk component. 
In the present paper we assume that in the very high state the innermost disk is deformed 
from an axisymmetric one to a two- armed one, 1 although the origin of the deformation is not 

1 In the warp resonant model (Kato 2004, 2008; Kato & Fukue 2006) the disk deformation was assumed to 
be one-armed. 
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explored. 2 In such deformed disks we impose a set of a two-armed (m = 2) vertical p-mode 
oscillation with negative energy and an axisymmetric (m = 0) g-mode oscillation with positive 
energy They have the same frequency. We then examine in what cases the set of the oscillations 
are resonantly excited through the disk deformation in order to apply the results to the HF 
QPOs in black-hole LMXBs. 

As is shown below, we consider two sets of the oscillations. The frequency ratio of 
these two set of oscillations derived by the present model is not always 3 : 2, depending on 
disk parameters. In the present model, magnetic fields in disks are important parameters to 
specify the frequencies of resonantly interacting oscillations. By adjusting the parameters in 
reasonable ranges, we examine whether we can describe the frequencies of twin HF QPOs 
of high-spin microquasars (GRS 1915+105, GRO J1655-40, XTE J1550-564, and H1743-322). 
One of important parameters is strength of poloidal magnetic fields, because the radial epicyclic 
oscillations are strongly modified by poloidal magnetic fields (Fu & Lai 2009). 

2. Outline of Present Model of HF twin QPOs 

2.1. Two-Armed-Deformed Disks 

As mentioned above, we assume that in the transition to the very high state the disk is 
deformed to a state with a two-armed pattern. The patern may rotate slowly in the azimuthal 
direction with frequency cjd- In this paper, for simplicity, the pattern is assumed just standing, 
= 0. (A generalization to a case of ud 7^ is simple.) We further assume that the two- armed 
deformation consists of two components concerning the node number(s) in the vertical direction. 
One is the deformation with one node in the vertical direction, i.e., = 1, and the other with 
no = 2. A riD = 1 deformation is the lowest mode of deformations which is asymmetric with 
respect to the equator. In summary, the sets of frequency, Wd, azimuthal wavenumber m D , and 
the vertical node number no of the deformation, (wd, ttld, wd), are both (0, 2, 1) and (0,2,2). 

2.2. Two Types of Oscillations and Their Propagation Regions 

Let us consider a set of oscillations. One is an axisymmetric (m = 0) g-mode oscillation, 
and the other is a two-armed (m = 2) vertical p-mode oscillation (see, e.g., Kato 2001 and 
Kato et al. 2008 for classification of disk oscillations). They are assumed to have the same 
frequency. These two oscillations can resonantly interact each other through the two-armed 
disk deformation, since the difference between the azimuthal wavenumbers is two and equal to 
mD- 3 This resonant interaction through the disk deformation can excite these two oscillations 

2 There are numerical simulations showing a deformation of the innermost part of disks. Machida & 
Matsumoto (2008) show that in a cool state of magnetized accretion disks, a torus is created in the in- 
nermost part of disks and deformed into a crescent-like shape. 

3 For the resonant interaction to occur, in addition to this condition concerning the difference of azimuthal 
wavenumbers, the difference of the vertical node numbers of the two oscillations must be 1 or 2. This is 
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simultaneously when wave energies of these oscillations have opposite sighs (Kato et al. 2011, 
Kato 2011b, see also Kato 2008). 

The axisymmetric g-mode oscillation has a positive energy because of the very fact of 
axisymmetry. The two-armed vertical p-mode oscillation has two propagation regions in the 
radial direction; inner and outer propagation regions. Between them, there are evanescent re- 
gion and the radius of corotation resonance. It is generally known that in the inner propagation 
region waves has a negative energy, while those in the outer propagation region have a posi- 
tive energy (e.g., Kato 2001). 4 Hence, if the propagation region of the axisymmetric g-mode 
oscillation and the inner propagation region of the two-armed vertical p-mode oscillation are 
spatially overlapped, we can expect simultaneous growth of these two oscillations with the same 
frequency by their resonant coupling through the disk deformation. The spatial overlapping of 
propagation region really occurs, as is shown below and in figures 1 and 2. 

Before discussing one more resonant condition which we adopt in this paper, we must 
notice the propagation regions of the two-armed vertical p-mode oscillation and of the axisym- 
metric g-mode oscillation, and their parameter dependences. First, we consider the m-armed 
vertical p-mode oscillations. At the boundary between the propagation and evanescent regions, 
the oscillations have long wavelength in the radial direction, and the relation between the fre- 
quency and boundary radius can be obtained by considering purely vertical oscillations (e.g., 
Kato 2011a). In the case of vertically isothermal disks with no magnetic fields, by solving the 
vertical eigen- value problem for purely vertical oscillations we find that the relation is given by 



(cf., Okazaki et al. 1987), where Q(r) is the angular velocity of disk rotation at the radius r, 
Q±(r) is the vertical epicyclic frequency, and n is a positive integer specifying the number of 
node of oscillations in the vertical direction, i.e., n = 1,2, ,... 5 

In the case of adiabatic oscillations in polytropic gases with polytropic index N, the 
radii where the vertical p-mode oscillations become purely vertical are 



[ &yfl 2 ± (n = 3) 

(Perez et al. 1997, Silbergleit et al. 2001, Kato 2005), where 7 = 1 + 1/JV. If the disk are 



related to the fact that tid is taken to tt-d = 1 and tid = 2. As will be discussed later, we consider axisymmetric 
g-mode oscillations with n = 1 and two-armed vertical p-mode oscillations with n = 2 and n = 3. Hence, this 
condition is satisfied. 

Let us consider an oscillation whose frequency is uj and the azimuthal wavenumber is m. Then, the corotation 
resonance occurs at the radius where uj — mfl. One of the propagaion region of the oscillation is inside the 
corotation radius, and the oscillation in the region has a negative energy since cj — mfl < there. 
In the case of vertical p-mode oscillations, n starts from 1. The case of n = is the p-mode. 




(1) 




(2) 
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vertically isothermal but terminated at a certain height, z s , by hot corona, and subject to 
toroidal magnetic fields which are distributed in the vertical direction so that c\/c 2 [c s (r) 
and ca(V) being respectively the acoustic and Alfven speeds due to the toroidal component 
of magnetic fields] is constant in the vertical direction, the radius where the vertical p-mode 
oscillations becomes purely vertical is determined by (Kato 2012b) 

2 2 

(^-^) 2 =(^f^^ + i)fii, (3) 

where K n ^ s is a value resulting from the eigenvalue problem of the purely vertical oscillations 
and depending on the node number, n, in the vertical direction and the dimensionless height 
of truncation, z s /H, H being the scale height of disks. The value of K n s + 1 is given in table 
1 of Kato (2012c). K n ^ s is n — 1 in the limit of no truncation and larger than n — 1 in general 
cases. 

As is shown above, the radii between the propagation and evanescent regions depends 
on the disk structure. Here, we write the boundary, for simplicity, in the form 

(u-mVL) 2 = i)tt\ (4) 

and if) is a parameter depending on disk structure. In the case of two-armed (m = 2) vertical 
p-mode oscillations, the boundary between the inner propagation and evanescent regions is thus 
given by 

u = 2tt-ip 1/2 {l ± . (5) 

In general, if) is larger than n [see equations (1) - (3)], and depends on n, c\/c 2 , and 
z s /H. In order to avoid complexity, we consider hereafter vertically isothermal disks with no 
truncation, i.e., z s = oo, and c\/c 2 is assumed to be radially constant. Then, if) is a function of 
b\ for a given n as [see equation (3)] 

^iTSP*"- 1 '* 1 - (6) 

where b 2 ^ is a parameter defined by = c\/c 2 . As is mentioned below, we are particularly 
interested in oscillations of n = 2 and n — 3. In the case of b\ = 0.4, we have if) = 2.17 for n = 2 
and if) = 3.33 for n = 3. Hence, in order to demonstrate typical examples of the frequency - 
radius relations given by equation (5) and (6), the relations in the case of if) = 2.17 and if) = 3.33 
are shown in figures 1 and 2 by thin curves for different values of a*. Figure 1 is for a* = 0.9 
and figure 2 is for a* = 0.5. The mass of the central source is taken to be 1OM in both figures. 
The left-hand side of the curves is the propagation region of the oscillations. 

Next, let us consider the propagation region of axisymmetric g-mode oscillations. The 
region in the propagation diagram (frequency - radius diagram) is given by u < k, when there 
is no poloidal magnetic fields in disks, and the boundaries of the propagation region of axisym- 
metric g-mode oscillation is specified by (Okazaki et al. 1987) 
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Radius (r/^) 

Fig. 1. Propagaton diagram (frequency - radius relation) for axisymmetric (m = 0) g-mode oscillations 
and two-armed (m = 2) vertical p-mode oscillations in the case of a* = 0.9 and M = 10M Q . Three thick 
curves give the relation of ui = k for three cases of b = 0, 0.1, and 0.2 from bottom to top. The propagation 
regions of axisymmetric g-mode oscillation in disks with 6 = 0, 0.1, and 0.2 are, respectively, below the 
corresponding curve. Among three thin curves the lower two curves are uj = 2Q — (2.17) 1 ^ 2 Q± (i.e., b\ = 0.4 
and n = 2) and uj — 2fl — (S^) 1 / 2 ^ (i.e., b\ = 0.4 and n = 3). Two-armed vertical p-mode oscillations 
with n = 2 and n = 3 can propagate, respectively, below the corresponding curve in the case of b\ = 0.4. 
The frequency given by the crossing point of two curves of uj = 251 — (2.17) 1 / 2 rij^ and uj = k is the frequency 
of the lower HF QPOs, ojl, in the present model. The waves with ojl propagate inside the point in the 
diagram (shown by arrow). Similarly, the frequency given by the cross point of uj = 2fl — (S^) 1 / 2 !!^ and 
uj = k is that of the upper HF QPOs. The waves with wh propagate inside the point in the diagram (shown 
by arrow). 




Fig. 2. The same as figure 1, except for a* = 0.5. The curves of uj = 2VL - (2.17) 1/2 ft± (i.e., fe| = 0.4 
and n = 2) do cross with the curves of u> = k in three cases of b = 0, 0.1, and 0.2. The curve of 
uj = 2Q — (3.33) 1 / 2 Q± (i.e., b\ = 0.4 and n = 3) crosses with the curve of u> = k in the case of b = 0. 
but not in other cases of b = 0.1 and 0.2. This shows the absence of ojh and HF QPO will be single. It is 
noted that in the case of b — 0, uj-r is smaller than cjl- Hence, in such cases the naming of wl and ojh due 
to the order of frequency should be changed. 
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UJ = K. 



(7) 



If poloidal magnetic fields exist, the radius specifying the boundary of propagation region 
of axisymmetric g-mode oscillations is strongly modified (Fu and Lai 2009). In the case of 
axisymmetric g-mode oscillations with n — 1, the propagation region is specified by u < k and 
k is given by (Fu and Lai 2009) 



where b = ca/c s , ca being the Alfven speed defined by the poloidal component of magnetic 
fields. The frequency - radius relation given by 



is also shown in figures 1 and 2 by thick curves. Three cases of b = 0, 0.1, and 0.2 are shown in 
each figure. The difference between figures 1 and 2 is the spin paramter, as mentiond before. 

2.3. A Working Hypothesis Determining Frequency of Resonant Oscillations 

The arguments given in subsection 2.2 are not enough to uniquely determine the fre- 
quency of resonantly excited oscillations, since the overlapping of the propagation regions is 
only a necessary condition. We should consider what case the growth of oscillations occur most 
strongly. To get a proper answer concerning this problem, detailed consideration on magni- 
tude of coupling terms will be necessary. Instead, we are satisfied here with adopting a simple 
working hypothesis, as was done in the case of the warp resonant model (Kato 2004, 2008). 

The excitation due to the resonance will be non-steady in the present problem in the 
sense that a wave will be excited at a place where the resonant interaction occurs most effi- 
ciently and then propagates in the radial direction. This wave will not satisfy in general the 
trapping condition. In other words, a possible condition determining the frequency of the ex- 
cited oscillations will be that at the same radius the group velocities of both oscillations vanish 
together so that the resonant interaction occurs efficiently. In the limit of waves of short wave- 
length, the group velocity vanishes at the boundary between the propagation and evanescent 
regions [see, for example, an expression for group velocity given by equation (12.16) by Kato 
et al. (2008)]. The boundaries are nothing but the radii of Lindblad resonances in the case of 
p- and g- mode oscillations. In the above contexts, we impose, as the condition determining 
the frequency of resonant oscillation, that the boundary radius between the propagation and 
evanescent regions is the same for two resonantly interacting oscillations. As mentioned before, 
the bounday between the propagation and evanescent regions is given by equation (5) for the 
two-armed vertical p-mode oscillations, and by equation (9) for the axisymmetric g-mode os- 
cillation. Thus, the crossing point of these two equations on the propagation diagram gives the 
frequency of the excited oscillations in the present model (see figures 1 and 2). 

In the followings, the above-mentioned condition is examined in two cases of a highly 
spinning source (a* = 0.9) and a less highly spinning source (a* = 0.5). First, we consider the 




(8) 



U) = hi 



(9) 
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coupling between a two-armed vertical p-mode oscillation with n = 2 and an axisymmetric g- 
mode oscillation with n = 1, taking a* = 0.9 (figure l). 6 The mass of the central source is taken 
to be 1OM . If if) = 2.17 (i.e., b\ = 0.4 and n — 2) and b = are adopted, the frequency is given by 
the crossing point of two curves of to = 2VL — (2.17) 1//2 S7j_ and to = k on the propagation diagram 
(see figure 1), the point being one of those labelled by cjl in figure 1. If b = 0.1 is adopted with 
other parameters unchanged, the point moves in the left-upper direction in the figure. If b = 0.2 
is adopted, the point moves further in the left-upper direction. The point is also labelled by 
wl in figure 1. Similarly, we consider the interaction between the two- armed vertical p-mode 
oscillation with n = 3 and the axisymmetric oscillation with n = l. 7 Now, we take ip = 3.33 (i.e., 
b\ = 0.4 and n — 3). Then, the crossing point of two curves of to = 2Vt — (3.33) 1/,2 f2j_ and to = k 
in three cases of b = 0, 0,1, and 0.2 is shown in figure 1. The points in two cases of b = and 
0.2 are shown by labelling wh. Here, we regards to-^ and uh in the same cases of parameters as 
the set of frequencies of the lower and upper HF QPOs. 

Next, we consider the case of a* = 0.5, keeping other parameters the same as those in 
corresponding cases of a* = 0.9. As is shown in figure 2, we have no crossing point between the 
curve of to = 20 — (3.33) ^Oj. (i.e., b\ = 0.4 and n = 3) and to = k in the region outside the 
radius of the marginally stable circular orbit, except for the case of b = 0. This comes from 
the following situations. Both curves of to = 2Vt — ip l l 2 VL^ and to = k shift downward on the 
propagation diagram as a* decreases, but the former much moves downward compared with 
the latter. It is noted that in the case of no poloidal magnetic fields, the set of wh) always 
exist, but cjl > wh- That is, for small a* beyond at a certain value, the role of ojl and wh is 
changed. 

We should notice that in the present model, HF QPOs do not always appear in pairs as 
mentioned above. That is, in the sources with high spin or weak poloidal magnetic fields, they 
can appear in pairs, but in other cases, there is the possibility of only one or no HF QPOs. 

3. Twin HF QPOs on Normalized Frequency — Frequency Diagram 

In order to compare the set of to^ and ton with observed frequencies of twin HF QPOs, 
we plot the set of (to-^M, to-nM) on the to^M - to^M diagram, where M is the mass of the central 
source. The reason why we consider the to^M - to^M diagram is that the position of calculated 
(wlM, to\{M) on the diagram is independent of M , and depends only on the parameters a*, b 
and b^- s For various sets of (a*, b\), the position of (to^M, to-^M) on the diagram is shown for 
three cases of b = 0, 0.1, and 0.2 (see subsequent subsections for individual source). Figure 3 is 
for b = 0.1, and figure 4 is for both 6 = and b = 0.2. The value of a* is changed continuously 

6 These two oscillations can resonantly interact through the disk deformation, since the difference of n's is 1 
and is equal to one of no- 

7 These two oscillations can resonantly interact since the difference of n's is 2 and is equal to one of nrj. 

8 If z s /H ^ oo, z s /H is also an additional parameter. 
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from a* = 0.99 towards a* = 0, while b\ is considered in four cases (b\ = 0, 0.2, 0.4, and 0.6) 
in figure 3 and in three cases (b\ = 0, 0.2, and 0.4) in figure 4. It is noted that the series of 
changing a* with b = 0.2 are terminated at certain value of a* (depending on b\) before goint 
to a* = (see the right-upper curves in figure 4). This is due to the absence of cjh (and wl too 
in some cases) as mentioned before (compare figures 1 and 2). 

Before evaluating (a*, b, b\) of four microquasars (GRS 1915+105, GRO J1655-40, XTE 
J1550-564, and H1743-322) by using the figures, we first draw a bird's eye view showing where 
twin HF QPOs of four microquasars are on this diagram. To do so, we adopt the following 
data for the microquasars: 

GRS 1915 + 105 : w£ bs = 113Hz, u^ s = 168Hz, M = 10 - 18M 

GRO J1655 - 40 : w° bs = 300Hz, u£ s = 450Hz, M = 5.1 - 5.7M 

XTE J1550 - 564 : w£ bs = 184Hz, u£ s = 276Hz, M = 8.5 - 9.5M 

H1743 - 322 : w£ bs = 166Hz, u^ s = 242Hz, M = 5.0 - 15.OM . (10) 

Observed lower and upper frequencies of the twin HF QPOs are denoted by a;£ bs and oj^ s in 
order to distinguish them from the calculated ones. Here, the observed mass ranges of GRS 
1915+105, GRO J1655-40, XTE J1550-564, and H1743-322 are taken, respectively, by referring 
to McClintock & Remillard (2006), Beer & Podsiadlowski (2002), Orosz et al. (2011), and 
Steiner et al. (2011), respectively. 

Because of uncertainty of mass, the position of (w£ bs M, o;g bs M) on the dia gram is on 
a finite line with a gradient of roughly 3:2 for each source. Since the ranges of the lines for 
four sources overlap on the diagram, we present the ranges by introducing rectangles whose 
diagonal shows the range. Using this convention, we overdraw in figures 3 and 4 the ranges of 
position of (w£ bs M, wg bs M) of the observed HF QPOs of the four microquasars by rectangles. 

It is noted that the position of twin QPOs by the epicyclic resonant model (Ep) (Kluzniak 
& Abramowicz 2001) and the warp resonant model (Wp) (Kato 2008) are also shown, for 
reference, in figure 3 for some values of a*. These models have no other parameter than a*. 
Hence, the position of twin HF QPOs by these models is determined only by a*. 

In figure 3 (and figure 4) the most probable points of (u;£ bs M , u^ s M) of the four mi- 
croquasars on the {u^M , u-rM) diagram are all roughly around the center of the figure. This 
suggests that the four microquasars belong all to a similar class. Furthermore, comparison of 
figure 3 with figure 4 shows that 6 = 0.1 is better than b = and b = 0.2 to describe observations. 
That is, figure 4 shows that if b = 0, the present model cannot describe observations even if 
a* is taken as high as a* ~ 1, unless each mass of microquasars is taken to be smaller than 
each mean value in the observationally allowed range given in equation (10). Figure 4 further 
shows that b = 0.2 is too large to describe observations in the same sense mentioned above. In 
other words, if the present model is correct, the poloidal magnetic fields in disks in these four 
microquasars will roughly be around b ~ 0.1. 
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Fig. 3. Relation between lo^M and ui^M on frequency - frequency diagram normalized by mass of the 
central source, i.e., the coordinates are frequency times mass. The relation is specified by three parameters, 
which arc a*, 6, and b\. The latter two are related to magnetic structure of disks, i.e., b is to the strength 
of poloidal magnetic fields and b\ to that of toroidal magnetic fields. In this figure b is taken to be b = 0.1. 
Four curves are for b\ = 0.0, 0.2, 0.4, and 0.6 (from left to right), a* being changed from 0.99 to lower 
values. The upper end of each curve is for a* = 0.99. Positions of four microquasars (GRS 1915+105, 
GRO J1655-40, XTE J1550-564, H1743-322) on the diagram are shown by rectangles. That is, due to 
uncertainties of mass, the position of (wlM, uj-rM) of each source on the diagram is on the line of diagonal 
of each rectangle specified by each source. For reference, the position of frequency - frequency relation by 
epicyclic model (Ep) and warp model (Wp) are also shown by open angles and filled angles, respecively. 
The straight line running from left-bottom to right-top is the line on which the frequency ratio is 3 : 2. 

Restricting only to the case of b = 0.1, we estimate the values of other parameters, 
especially the spin parameter a* in order to examine whether it is consistent with that estimated 
from spectral analyses. Hereafter, we discuss the four microquasars separately. 

3.1. GRS 1915+105 

Figure 5 is an extension of a part of figure 3. The diagonal of the large rectangle covering 
the alomst all area of figure 5 shows the range of the possible position of the set of observed 
twin QPOs, when the mass range is taken to be M = 10 — 18M Q . The small filled rectangles on 
the diagonal in figure 5 is the position of the twin QPOs in the cases where the mass is taken 
from 1OM to 18M & . The four lines running in the middle region of the figure are for b\ = 0.0, 
0.2, 0.4, and 0.6 from left to right, changing the value of a*. On each line the upper end point 
is for a* = 0.99. This figure shows that the set of (a*, b\) which can describe the observed 
set of QPO frequencies are (~ 0.93, ~ 0.3) if M = 14M , and (0.99, - 0.7) if M = 16M . 
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Fig. 4. The same as figure 3 except that in this figure the frequency - frequency relation in two cases of 
b = and b = 0.2 are shown. 




Fig. 5. An extension of figure 3 so that we can inspect what set of (a*, b\) well describe the observed 
twin HF QPOs of GRS 1915+105 in the case where b = 0.1 is adopted. 

If M > 16M , the present model cannot describe the observed twin QPOs as long as b = 0.1 
is adopted, since the value of a* required becomes larger than unity. In other words, if this 
source has really M > 16M and the present model is correct, the poloidal magnetic fields of 
this source are slightly stronger than b = 0.1. It is noticed that if b is larger than b = 0.1, the 
curves on figure 5 shift in the right-upward direction, and observations can be described by a* 
smaller than unity (compare figures 3 and 4). A large value of a* presented here is consistent 
with a* > 0.98 derived from spectral analyses (McClintock & Remillard 2006). 
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Fig. 6. An extension of figure 3 so that we can inspect what set of (a*, b\) well describe the observed 
twin HF QPOs of GRO J1655-40 in the case where b = 0.1 is adopted. 

3.2. GRO J1655-40 

Figure 6 is for GRO J1655-40, and the same as figure 5 except that the twin QPO 
frequencies are taken here to be w£ bs = 300Hz and w* s = 450Hz. Figure 6 shows that the set of 
(a*, b\) which can describe observations are (~0.90, 0.2) when M = 5.1M , and (~0.95, ~0.4) 
when M = 5.5M Q , and (~ 0.99, ~ 0.5) when M = 5.7M . The spin estimated from spectral 
analyses is a* = 0.70 - 0.80 (Shafeee et al. 2006) and a* ~ 0.98 (Miller et al. 2011). Our results 
are closer to the latter, as long as b = 0.1 is adopted. If b is slightly larger than b = 0.1, however, 
a smaller a* (with a smaller b%) fits observations (compare figures 3 and 4). 

3.3. XTE J1550-564 

Figure 7 is for XTE J 1550-564, since w£ bs = 184Hz and wg bs = 276 Hz are adopted. 
This figure shows that the set of (a*, b\) which can describe observed twin QPO frequencies 
are (~0.92, ~ 0.2) when M = 8.5M , and (~ 0.96, 0.4) for M = 9.1M Q . The value of spin 
parameter a* required by the present model is rather high compared with the value derived 
from spectral analyses. The latter is a* ~ 0.5 in average (Steiner et al. 2011) or a* = 0.70 — 0.75 
(Miller et al. 2009). 

This difference will be accounted for if this source has a stronger poloidal magnetic fields 
than GRS 1915+105 in the sense that b > 0.1. This is because if b > 0.1, the calculated curves 
in figure 7 shift in the upper-right direction in the figure, as mentioned before in relation to 
GRO J1655- 40 (compare figures 3 and 4). 

3.4. H1743-322 

Figure 8 is for H1743, since this diagram is drawn by taking o;£ bs = 166 Hz and =242 
Hz. This figure shows that the set of (a*, b%) which can describe observations are (~0.80, ~ 0) 
for M = 8M , (~0.9, ~0.2) for M = 9M G , and (0.93, 0.4) for M = 1OM . The present model 
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Fig. 7. An extension of figure 3 so that we can inspect what set of (a*, b\) well describe the observed 
twin HF QPOs of XTE J1550-564 in the case where b = 0.1 is adopted. 
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Fig. 8. An extension of figure 3 so that we can inspect what set of (a*, b\) well describe the observed 
twin HF QPOs of H1743- 322 in the case where b = 0.1 is adopted. 

cannot describe observations if M > 11M Q , as long as b = 0.1 is adopted. Steiner et al. (2011) 
estimate a* to be smaller than 0.7 by spectral analyses. This difference of our a* and that of 
Steiner et al. (2011) can be accounted for if stronger poloidal magnetic fields than b = 0.1 are 
present in H1743, as in the case of XTE J1550-564. 

4. Discussion 

A standpoint of this model is that the 3 : 2 frequency ratio observed in twin HF QPOs 
of microquasars does not come from their excitation processes, but occurs by chance. In the 
present model, a quasi-periodic oscillation is a result of a resonant interaction between two 
oscillations with positive and negative energies through disk deformation (Kato et al. 2011). 
As the disk deformation, a two-armed (mo = 2) deformation is assumed, different from our 
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previous studies where an one-armed (rriD = 1) deformation is considered (e.g., Kato 2004, 
2008). As the set of positive and negative energy waves, we consider an axisymmetric (m = 0) 
g-mode oscillation with positive energy and a two-armed (m = 2) vertical p-mode oscillation 
with negative energy. 

The difference of azimuthal wavenumber in the above two disk oscillations is two, which is 
equal to the azimuthal wavenumber, mo, of the disk deformation. Hence, a necessary resonant 
condition among azimuthal wavenumbers of the two oscillations and the disk deformation is 
satisfied. For the resonance really to occur, a relation among vertical node number is also 
necessary for the two oscillations and the disk deformation. In this paper we have considered 
two two-armed vertical p-mode oscillations with two nodes (n = 2) and three nodes (n = 3). On 
the other hand, as an axisymmetric g-mode oscillation, we have considered the mode with one 
node (n = 1) in the vertical direction. Hence, for these vertical p-mode oscillations to resonantly 
interact with the g-mode oscillations through the disk deformation, the disk deformation must 
have both components of n — 1 and n = 2. The deformation with n — 1 is the fundamental 
deformation of disks which is anti-symmetric with respect to the equator. 

One of vague points of this model is how much the assumption concerning the frequency 
of excited oscillations (i.e., coincidence of boundaries of propagation regions for two interacting 
oscillations) is relevant or not. To study this problem in details, we must examine functional 
forms in the radial and vertical directions of two interacting positive and negative energy 
oscillations and evaluate the volume integration of coupling terms by using those functional 
forms of oscillations and disk deformation [see Kato (2008) as an example of detailed calculations 
of the coupling terms]. This is a complicated problem and beyond the scope of this paper. 

Different from the kHz QPOs in neutron-star LMXBs, the observed frequencies of twin 
QPOs of black-hole LMBXs are robust, they being time-independent and their ratio being close 
to 3 : 2. Furthermore, for four microquasars considered in this paper, the observational points 
of (u)ft S M, cu£ bs M) on the (curM, u^M) diagram are close when the typical observational values 
of M's are adopted for each source (see that the centers of four rectangles in figure 1 are close). 
These observational evidences suggest that the disks of these sources are roughly in a similar 
state with similar time-independent values of disk parameters. In the present model we have 
two dimensionless disk parameters describing the magnetic structure of the innermost region 
of disks. They are b and b^, where b = ca/c s and by = ca/c s , ca and ca being the Alfven 
speed defined by the poloidal and toroidal magnetic fields in disks, respectively, and c s being 
the acoustic speed. As discussed in the text, b > 0.1 and b T = (0.2) 1 / 2 ~ (OA) 1 / 2 = 0.45 ~ 0.63 
seem to be relevant to describe observations. Detailed comparison of the present model with 
observations will not be instructive at the present stage, since the validity of the working 
hypothesis mentioned in subsection 2.3 might not be clear enough. The strength of magnetic 
fields required to describe the QPO frequencies of the four microquasars by the present model is, 
however, in a narrow range as noticed above. This may suggests that the strength of magnetic 
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fields in disks is self- adjusted by a balance between winding by differential rotation and release 
by jet ejection. 

In the present model, QPOs are not always in pairs. In some sources there is only one 
QPO or no QPO, depending on parameters. For example, in disks with b = 0.2 and b\ = 0.4, 
two QPOs (n — 2 and n — 3) are present when a* = 0.9 (see figure 1), but QPO is only one 
(n = 2) if a* is as small as 0.5 (see figure 2), because the curve of u = 2Vt — ip 1 / 2 ^}^ for n = 3 
(with b\ = 0.4) is below the curve of u = R (with b = 0.2) in all radius. Even the QPO of n = 2 
disappears (there is no QPO) if b\ is slightly larger than b\ = 0.4 {b = 0.2 and a* = 0.5), as 
understand from figure 2. Next, let us consider sources with b — 0. In this case, there are always 
two QPOs as shown in figures 1 and 2. However, in the case of figure 2 (a* = 0.5), the cross point 
between tu = k and u = 2VI - ^/ 2 fi ± occurs inside the radius where k becomes the maximum, 
as is shown in figure 2. In this case the propagation regions of the two oscillations are in the 
oposite side of the crossing point, and there is no overlapping region of the two oscillations. In 
such cases, the coupling of the two oscillations and thus excitation of the oscillations will be 
weak. In summary, we can say that the sources where twin QPOs are observed are those where 
the disk parameters are in some limited range with strong spin parameter. This might be one 
of reasons why twin QPOs are observed only in limitted sources. 

Finally, we should note that we have an important disk parameter which was not con- 
sidered in this paper, but will affect results. This is the vertical thickness of disks, z s . In this 
paper we assumed that the disk is isothermal in the vertical direction and extends infinitely. In 
actual situations, the disks will be truncated in the vertical thickness at a certain height, z B , by 
the presence of corona. In such disks, the propagation region of the two-armed vertical p-mode 
oscillations is modified (Kato 2012b, c), which changes wh and wl- This will be a problem to be 
examined. It is noticed that retreat of the inner edge of disks from the radius of the marginally 
stable circular orbit has no serious effects on results of the present model, since the oscillations 
treated in the model are not standing ones and no inner boundary condition is imposed to the 
oscillations. 

A note added on July 24, 2012. 

In order to examine validity of the wave excitation mechanism adopted in this paper 
[i.e., excitation of a set of waves with positive and negative wave energies through a resonant 
coupling with disk deformation (Kato et al. 2011)], we consider waves in tidally deformed disks 
of binary stars. It is found that the tidal instability and superhump phenomena in dwarf novae 
can well be interpreted by the mechanism. This will be shown in a subsequent paper. 
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